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Abstract 



^ I In the paper ( math-ph /0504049| Jarlskog gave an interesting simple parametriza- 

tion to unitary matrices, which was essentially the canonical coordinate of the second 
kind in the Lie group theory (|math-ph / 0505047 ) . 

In this paper we apply the method to a quantum computation based on multi- 
level system (qudit theory) . Namely, by considering that the parametrization gives a 
complete set of modules in qudit theory, we construct the generalized Pauli matrices 
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which play a central role in the theory and also make a comment on the exchange 
gate of two~qudit systems. 

Moreover we give an explicit construction to the generalized Walsh-Hadamard 
matrix in the case of n = 3, 4 and 5. For the case of n = 5 its calculation is relatively 
complicated. In general, a calculation to construct it tends to become more and 
more complicated as n becomes large. 

To perform a quantum computation the generalized Walsh-Hadamard matrix 
must be constructed in a quick and clean manner. From our construction it may be 
possible to say that a qudit theory with n > 5 is not realistic. 

This paper is an introduction towards Quantum Engineering. 

1 Introduction 

In the paper 1] Jarlskog gave a recursive parametrization to unitary matrices. See also 
[2] as a similar parametrization. One of the authors showed that the recursive method 
was essentially obtained by the so-called canonical coordinate of the second kind in the 
Lie group theory, jS]. 

We are working in Quantum Computation, therefore we are interested in some appli- 
cation to quantum computation. 

One of key points of quantum computation is to construct some unitary matrices 
(quantum logic gates) in an efficient manner like Discrete Fourier Transform (DFT) when 
n is large enough, j^. However, such a quick construction is in general not easy, see |3] 

or i, m 

The parametrization of unitary matrices given by Jarlskog may be convenient for our 
real purpose. We want to apply the method to quantum computation based on multi-level 
system (qudit theory). One of reasons to study qudit theory comes from a deep problem 
on decoherence (we don't repeat the reason here). In the following let us consider an n 
level system (for example, an atom has n energy levels). 

Concerning an explicit construction of quantum logic gates in qudit theory, see for 
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example [S], jH] and JU]- By use of the new parametrization to unitary matrices we want 
to construct important logic gates in an explicit manner^, especially the generalized 
Pauli matrices and Walsh-Hadamard matrix, which play a central role in qudit theory. 

In this paper we construct the generalized Pauli matrices in a complete manner, while 
the Walsh-Hadamard matrix is constructed only for the case of 3, 4 and 5 level systems. 
A calculation to construct it for the 5 level system is relatively complicated compared 
to the 3 and 4 level systems. In general, a calculation tends to become more and more 
complicated as n becomes large. 

The generalized Walsh-Hadamard matrix gives a superposition of states in qudit the- 
ory, which is the heart of quantum computation. It is natural for us to request a quick 
and clean construction to it. 

Therefore our calculation (or construction) may imply that a qudit theory with n>b 
is not reahstic. Further study will be required. 

2 Jarlskog's Parametrization 

Let us make a brief introduction to the parametrization of unitary matrices by Jarlskog 
with the method developed in 0. The unitary group is defined as 




(1) 



and its (unitary) algebra is given by 




(2) 



Then the exponential map is 




(3) 



This map is canonical but not easy to calculate. 



Quantum computation is not a pure mathematics, so we need an explicit construction 
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We write down the element X e u{n) explicitly : 



X = 





Zl2 


Zl3 
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ie 



This X is decomposed into 



X = Xo + X2 + --- + X. + --- + Xn 



where 



ien-i 



ie, 



and for 2 < J < n 



-{Zi 



\ 







Z2j 



(4) 



(5) 



(6) 



Then the canonical coordinate of the second kind in the unitary group (Lie group) is 
well-known and given by 



u{n) 3 X ^ Xq + X2^ h H h X„ 
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e^^ . . . e^" e C/(n) (7) 



in this case ^. Therefore we have only to calculate e^^ for j > 2 (j = is trivial), which 
is easy (see Appendix). The result is 

( - (l - cos{^ {zj\zj))) \zj){zj\ {zj\zj))\zj) ^ 



e^^ 



-sm{^{zj\zj)){zj\ 



C0S{^ {Zj\Zj)) 



\ 



where \zj) is a normalized vector defined by 
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= — ki) =^ = 1- 



(9) 



We make a comment on the case of n = 2. Since 



|5) 



e^°, (51 



= = 1, 



we have 



e^°e^2 
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cos(|z|) e*"sin(|z|) 



— e *° sin(|2;|) cos(|2;|) 



/ gjQ/2 
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cos(|z|) sind^l) \ I e-*"/2 



^ g-W2 j y -sin(|z|) cos(|z|) 



cosd^l) sin(|2;|) 



\ 



Ji92-a/2) 



Ja/2 



(10) 



-sin(|2;|) cos(|2;|) / y 

This is just the Euler angle parametrization. 

Therefore the parametrization (|7j) may be considered as a kind of generalization of 
Euler 's angle one. 



In the following we set 



Ao = e 



Xo 



V 



J,dn 



^ There are of course some variations 



(l 



- (1 - COS Pj) 1%) I sin 
— sin Pj{zj\ cos Pj 




(11) 



V 




for 2 < j < n. More precisely, we write 



Ao = Ao{{9i,92,- ■ ■ ,0n}), Aj = Aj{{zij,Z2j,- ■ ■ ,Zj-ij};Pj) for j = 2,---,n (12) 



including parameters which we can manipulate freely. 

From now on we consider Aj a kind of module of qudit theory for j = 0, 2, ■ ■ ■ , n, 
namely {Aj\ j = 0,2, ■ ■ ■ ,n} becomes a complete set of modules. By combining them 
many times ^ we construct important matrices in qudit theory in an explicit manner. 



Let us make a brief introduction to a qudit theory. The theory is based on an atom with 
n energy levels {(|/c), Ek) | < A; < n — 1} , see the figure 1. 

First of all we summarize some properties of the Pauli matrices and Walsh-Hadamard 
matrix, and next state corresponding ones of the generalized Pauli matrices and general- 
ized Walsh-Hadamard matrix within our necessity. 

Let (72, cTs} be Pauli matrices : 



By (fTT^ (72 = io'io'3, so that the essential elements of Pauli matrices are {ai, a^} and they 
satisfy 



3 Qudit Theory 




(13) 



2 2 T t t ITT 



(14) 



A Walsh-Hadamard matrix is defined by 




(15) 



•^we take no account of an ordering or a uniqueness of {Aj} in the expression Q 
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Figure 1: Atom with n energy levels 



This matrix is unitary and it plays a very important role in Quantum Computation. 
Moreover it is easy to realize it in Quantum Optics as shown in jH]. Let us list some 
important properties of W : 

W'^ = I2, = W = W~\ (16) 
ai = Wa3W-\ (17) 

The proof is very easy. 

Next let us generalize the Pauli matrices to higher dimensional cases. Let {Ei, S3} be 
the following matrices in M(n, C) 



Si 





1 
1 

1 • 



\ 



1 



1 



a 



a 



J V 
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where o" is a primitive root of unity cr"^ = 1 (cr = e n ). We note that 



a = (7"-\ l + (j+--- + (j"-^-0. 

Two matrices {Si, E3} are generahzations of the Pauh matrices {ai, a^}, but they are not 
hermitian. Here we hst some of their important properties : 



^1 — ^3 — -^n , ^1—^1 ) ^3 — ^3 ) ^3^1 — Cr^l^3 ■ 

For n — 3 and n — 4 Tii and its powers are given respectively as 



Si 



1 
1 



y2 
> ^1 
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1 



and 
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1 
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^0 1 ^ 
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y3 



1 

1 
1 

1 



If we define a Vandermonde matrix based on a as 

^ 1 1 

1 



1 ^2 
1 (7 



1 

r(n-l)(n-2) 



r2(n-l) 



(19) 



(20) 



(21) 



(22) 
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ll 1 1 




1 \ 


1 a 






1 




0-2(n-l) 


1 cr"-2 o-^(""^) 




^{n-l)(n-2) 


1 n— 1 2(n— 1") 




' J 


w'^w = WW^ = 


In, 




Si = wj:3W^ = 




1 



(23) 



then it is not difficult to see 

(24) 
(25) 

Since W corresponds to the Walsh-Hadamard matrix (fTHjl . it may be possible to call 
W the generalized Walsh-Hadamard matrix. If we write W"^ = (wah), then 

1 1 /27ri \ 

Wab = —;=o'"'^ = ^^exp ( ab ] for < a, b < n — 1. 

i/n \ n } 

This is just the coefficient matrix of Discrete Fourier Transform (DFT) if n = 2^^ for some 
A; G N, see 

For ra = 3 and n = 4 is given respectively as 



73 



^ 1 1 ^ 



1 a' a 
1 a 



/ 



1 
71 
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(26) 



and 



/i 1 1 1 ^ 



,3 ^2 

1 1 



1111 

1 -Z -1 i 

1-11-1 
1 Z -1 -Z 



(27) 



We note that the generalized Pauli and Walsh-Hadamard matrices in three and four 
level systems can be constructed in a quantum optical manner (by using Rabi oscillations 
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of several types), see and 0. Concerning an interesting application of the generalized 
Walsli-Hadamard one in three and four level systems to algebraic equation, see [TT] . 



4 Explicit Construction of the Generalized Pauli and 
Generalized Walsh— Hadamard Matrices 

First let us construct the generalized Pauli matrices. From (fT^ it is easy to see 

Ao({0, 27r/n, 47r/n, ■ ■ ■ , 2(n - l)7r/n}) = S3. (28) 



Next we construct Si. From (fT2ll we also set 



A,- = /l,({0,---,0,l};7r/2) 



for j = 2, ■ ■ ■ ,n. Then it is not difficult to see 

' 



A2A3 ■■■An 



1 
-1 



-1 
0-10 



■1 



(29) 



(30) 



Therefore if we choose Aq as 



= Ao({0,7r,---,7r}) 



-1 



-1 



(31) 



then we finally obtain 



Aov42A3---A„ = Si. 



(32) 
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I a) 



\b) 



\b) 



I a) 



Figure 2: Exchange gate on two-qudit system 



From ()28p and ()32|) we can construct all the generalized Pauli matrices 

|gj0{a,fe)^a^fe |0<a, & — l|, where e*"^*-"'*-* is a some phase depending on a and 6. 

Similarly we can construct the matrix 



K 



(33) 



as follows. If n = 2k, then 

Ao({0^_^, 7V;^})Afc+2({0, ■ • ■ , 0, 1, 0}; 7r/2)A,+3({0, • ■ ■ , 0, 1, 0, 0}; 7r/2) ■ • • 
k+1 k-1 

xA2k-i{{0, 0, 1, 0, ■ ■ ■ , 0}; 7r/2)A2fc({0, 1, 0, ■ ■ ■ , 0}; n/2) = K (34) 
and if n = 2A; — 1, then 



^({0^_^, 7v;^})^fe+i({0, ■ ■ ■ , 0, 1}; 7r/2)Afc+2({0, ■ ■ ■ , 0, 1, 0}; 7r/2) • • • 
k k-1 

X A2fc_2({0, 0, 1, 0, ■ ■ ■ , 0}; 7r/2)A2fc„i({0, 1, 0, ■ ■ ■ , 0}; n/2) = K. (35) 

Both Si and K play an important role in constructing the exchange (swap) gate in 
two-qudit systems like the figure 2 where S=Si. To be more precise, see |12j . 
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It is interesting to note the simple relation 



= K. (36) 
Namely, the generalized Walsh-Hadamard matrix W (j^^ is a square root of K. 



Second we want to construct the generalized Walsh-Hadamard matrix, which is how- 
ever very hard. Let us show only the case of n = 3, 4 and 5. 



(a) = 3 : For (j^UI) we have 



where each of matrices is given by 



W = A0A3A2AQ 



v4o = v4o({0, 277/3, 47r/3}) 



J2tt/3 



( g-i7r/12 



a; = A;({-vr/12, 777/12,0}) 



-,i47r/3 



^j77r/12 



and 



A3 = A3 ({I/V2, 1/V2}; cos-i(l/V3)) = 



^/3+l 


v^-1 


1 


\ 


2 


2 




\/3-l 


v/3+1 


1 




2 


2 




-1 


-1 


1 


/ 



and 



A2 = A2 ({e-*-/2};7r/4 



i 1 
V2 V2 



Here we have used 



cos(7r/12) = and sm(7r/l/j = 



4 



In this case, the number of modules is 4. 



(37) 
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(b) n = 4 : For ()27|) we have 



W = AoA4SA3A2A'qS 



where each of matrices is given by 



Ao = Ao({0,27r/4,47r/4,67r/4}) 



A; = ylo({7r/4,57r/4,0,0}) 



V2 



1+1 
V2 



and 



and 



^4 = ^4 ({l/VS, l/VS, l/VS}; 7r/3 



A3 = A3 ({1/72, 1/V2}; cos"i(-l/3) 



/ 5 

6 

_ 1 
6 

_ 1 

6 

1 

V -2 



/ 1 

3 

_2 
3 

_2 

3 



and 



5 = Ao({0,0,7r,0})A3({0,l};7r/2) 
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and 



A2 = A2 ({e'-/2|.^/4 



a/2 ^/2 

i 1 
V2 V2 



\ 



J 



In this case, the number of modules is 9. 

Last we show a calculation for the case of n = 5. However, it is relatively complicated 
as shown in the following. 



(c) n = 5 : We have 



W = AoA,A^SiAsS2A'^ 



where each of matrices is given by 



(39) 



Ao = Ao{{0, 2n/5, 47r/5, Gtt/S, Stt/S}) = 



= ylo({97r/10, IStt/IO, -Stt/IO, tt/IO, 0}) 
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where 



and 



a — e 



i27r/5 



cos(27r/5) + isin(27r/5) 



V5-1 .V10 + 2V5 



+ ^ 



^ = A, ({1/2, 1/2, 1/2, 1/2}; cos-i(l/V5 
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. /.X . / /.N Vl0-2\/5 .\/5 
= sin(27r/5), a = sin(7r/5) + z-— = -^^ cos(^4) 

^ 4r ^ 

r 2(35+^)/, sin(2ir/5)\ 



Si = A„({0, w, TT, 0, 0})AM1}; w/2)As({0, 1}; 7r/2) = 



^ 1 ^ 



1 
1 



52 = ^o({0,0,7r,0,0})/l3({l,0};7r/2) 



^ 1 ^ 



1 
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and 



^3 = A3 ({-l3/v,'P/v}-9^ 



V^(v^-(a+ta2)) v^(y5-(a+ta2)) 

1 _ W A. 

V5(V5-(a+ta2)) V5( V5-(a+ta2)) ^ 



V5 



V5 



V5 



g+tffl 
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where 



/? = a + taa, = y 5 — (a + ta^)^ = \/2|/3|, cos( 
^ ^ 2(7x/5 + l) / ^ sin(27r/5)\ 



a + ta 
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A comment is in order. 

(1) Our construction is not necessarily minimal, namely a number of modules can be 
reduced, see [12]. 

(2) For n > 6 we have not succeeded in obtaining the formula like (jHTjl or or In 
general, a calculation tends to become more and more complicated as n becomes large. 

(3) The heart of quantum computation is a superposition of (possible) states. Therefore 
a superposition like 

|o)- 



|0) + |1) + |2) + |3) + |4) 



in the 5 level system must be constructed in a quick and clean manner. The generalized 
Walsh-Hadamard matrix just gives such a superposition. Our construction in the system 
seems to be complicated, from which one may be able to conclude that a qudit theory 
with n > 5 is not realistic. 



5 Discussion 

In this paper we treated the Jarlskog's parametrization of unitary matrices as a complete 
set of modules in qudit theory and constructed the generalized Pauli matrices in the 
general case and generalized Walsh-Hadamard matrix in the case of n = 3, 4 and 5. 
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In spite of every efforts we could not construct the Walsh-Hadamard matrix in the 
general case, so its construction is left as a future task. 

However, our view is negative to this problem. In general, a calculation to construct 
it tends to become more and more complicated as n becomes large. See ()37|1 and ()H8|1 and 
fl39|) . Therefore it may be possible to say that a qudit theory with n > 5 is not realistic. 
Further study will be required. 

Our next task is to realize these modules {Aj\ j = 0, 2, ■ ■ ■ , n} in a quantum optical 
method, which will be discussed in a forthcoming paper. 
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Appendix Proof of the Formula 



In this appendix we derive the formula (jH)) to make the paper self-contained. 
Since 



■(^1 



V 



(k 

V On . 



e^^ 



we have only to calculate the term , which is an easy task. From 



K 



-{z,\ 



^2 I 



-{Zj\Zj) ) 



I 



Oj-1 -(^iki)ki) 



we have important relations 



(40) 



i^2-+^ = (-(z,|^,))"K, K'-^^ = {-i,z,\z,)rK^ for n>0. 
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Therefore 



OO 1 OO -| 



00 -j^ 



00 1 



n=0 



(2n + 1)! 



n=0 



(2n + 2)! 



n=0 
2n+2 



(2n + 1)! 

1 



, 2n+l 



(2n+ 1)! 



, 2n+l 



1, - 



1 



(2n)! 



,/{^to (2n+l)! 



Ij + ( 1 — cos 



If we define a normalized vector as 



then 



{Zj\Zj) 



I -(5,1 



{Zj\zj) 



Therefore 



Ij-i - (1 - cos(y(I^)) sm{^ {zj\zj))\zj) 

C0S{^ {Zj\Zj)) 



As a result we obtain the formula (jH)) 

^ Ij-i - (1 -cos(^(^j|^j))) \zj){zj\ sm{^ {zj\zj))\zj) ^ 
-sm{^ {zj\zj)){Sj\ cos{^ {zj\zj)) 



e^^ 



(41) 



(42) 
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